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Abstract 

We describe the application of renormalization group improved perturbative QCD 
to inelastic lepton-hadron scattering at high center-of-mass energy but compar- 
atively low photon virtuality. We construct a high energy factorization theorem 
which complements the mass factorization theorem used for processes with high 
virtualities. From it we derive a renormalization group equation which resums all 
large logarithms at high energy, thereby extending to this regime asymptotic free- 
dom and thus the full range of perturbative computational techniques. We discuss 
the solution of this equation in various limits, and in particular show that the high 
energy behaviour of physical cross-sections is consistent with phenomenological 
expectations and unitarity bounds. 



March 1997 

* Royal Society University Research Fellow 



Asymptotic freedom Jl]] in deep inelastic processes, derived using mass factorization 
and the renormalization group, lies at the heart of our understanding of perturbative 
QCD. Attempts to use perturbative QCD to understand processes with high centre-of- 
mass energies but relatively low virtualities are as old but have never had as firm a 
theoretical foundation. In these processes distances from the light-cone are no longer short, 
but longitudinal distances are compressed by the large relative longitudinal momenta of 
the interacting particles ||. At leading order in perturbation theory cross-sections are 
asymptotically flat, while at higher orders logarithms of energy (or rapidity) generate a 
rise. However these large logarithms eventually spoil the perturbative expansion. Attempts 
to sum them |||| lead to cross-sections which seem to rise too quickly to agree with current 
data [0, and eventually violate unitarity bounds |J. 

Here we describe an approach to these processes based on a high energy factorization 
theorem for inelastic lepton-hadron scattering which, when combined with the renormal- 
ization group, implies that QCD becomes asymptotically free at high energies. The large 
energy logarithms can then be resummed to give a self-consistent perturbative description 
of inelastic scattering at high energy. A partonic picture is developed similar in many 
ways to the conventional one useful at high virtualities. We find that cross-sections 
then rise only logarithmically, consistent both with phenomenological expectations and 
the behaviour of unitarity bounds. In this paper we will outline our main results; detailed 
proofs and a wider range of applications will be presented in a forthcoming publication. 

Perturbative logarithmic corrections to cross-sections arise (in physical gauges) from 
two particle reducible diagrams: two particle irreducible diagrams can be at most constant 
at high scales in a nonabelian gauge theory because they contain no logarithmic singular- 



ities [HJ. The two particle reducible contributions to an inelastic virtual photon-hadron 



cross-section satisfy the fundamental factorization [DJ 



1 2 

a^,QV^) = [^£ d £c(l,^a s ^))f(y,k*-,S), (1) 

where, if q and p are the momenta of the virtual photon and hadron, Q 2 = —q 2 , x = 
Q 2 /2p ■ q, A is a fixed scale, \i is the renormalization scale, and parton indices have been 
suppressed. This factorization may be proven formally in axial gauges by a projection 
argument [nj which generalises that used to prove mass factorization [1(],|12| . The content 



of eq. (|l|) is that all collinear singularities can be factored into the unintegrated two particle 
irreducible parton distribution function / in such a way that the process-independent two 
particle reducible coefficient function C only depends on the renormalization scale through 
the strong coupling. The dependence on the parton and hadron momenta is then fixed by 
simple kinematics. Here we will take ([!]) as our starting point, show how to derive mass 
factorization from it and then, by a parallel argument, derive the energy factorization 
theorem which will be the basis for our high energy perturbative expansion. 

To expose the symmetry implicit in ([]]), we introduce S 2 = A 2 /x, I 2 = A 2 /?/, so that 

{2) (S 2 q 2 2 \ rdi 2 rdk 2 ( s 2 q 2 2 \ / 1 2 p A 
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where the full dependence on the renormalization scale /i 2 is displayed explicitly (though 
that of the cross-section is of course fictitious). The formal extension of the integration 
limits is of no consequence for the applications we shall consider (provided the definitions 
of coefficient functions and parton distributions are suitably extended). Just as the high 
virtuality limit corresponds to Q 2 — > oo at fixed S 2 , so the high energy limit is now S 2 — ► oo 
at fixed Q 2 . Both the convolutions may be undone by taking Mellin transforms: defining 
for example 

(2) f 2s [°°dS 2 f A 2 \N f°° dQ 2 /A 2 \ M ( S 2 Q 2 A 
°nm(p)= -02-02 / -7^-7^2 I 772"' TT2"' ^ ) ' ( 3 ) 



J S 2 \S 2 J J Q 2 \Q 2 J V ^2 
the double factorization @ assumes the simple form 

a NM = C NM (a s (p 2 ))f NM (v> 2 )- (4) 



We next explain how mass factorization can be derived from the double factoriza- 
tion (|^). At large Q 2 we may write ([2]) as 

a%\Q 2 /v 2 ;fi 2 ) = — dM(?L) C NM [a 8 {n 2 ))f NM {n 2 ), (5) 

with the integration contour passing to the right of the origin and closed on the left. 
The singularities in the M-plane at any fixed order in perturbation theory may be 
found by power-counting arguments [[13|, since the Mellin transform of (A 2 /Q 2 ) k ln n is 



n!/(M — k) n+1 . It follows that the leading singularity at large Q 2 is at M = 0, singular- 
ities further to the left giving rise to higher twist terms which vanish as powers of 1/Q 2 - 
Similarly all two particle irreducible contributions to the cross-section are down by powers 
of l/Q 2 , and can thus be ignored at leading twist. 

Hence at large Q 2 and leading twist we are interested only in the behaviour of Cnm 
and fNM in the neighbourhood of M — 0. The two particle irreducible unintegrated parton 
distribution function is regular at M = (again by power counting ||13|| ), while the (two 
particle reducible) coefficient function contains logarithmic singularities which correspond 
to multiple poles at M = 0. We may thus perform a Laurent expansion of the coefficient 
function 

00 

C NM (a a (jj?)) = C%{a 8 {i?))M- m -\ (6) 

m= — 00 

so that at large Q 2 

00 

C N (Q 2 /k 2 -^ 2 ) = -r^(« s (^ 2 ))(lnQ 2 AT + 0(k 2 /Q 2 ). (7) 

Since the logarithmic singularities are generated perturbatively, the coefficient C^(a s ) will 
begin at 0{a™) when m > 0. Similarly the unintegrated parton distribution function has 
a Taylor expansion 

00 

fNM^ 2 ) = Yl fN(^)M m , (8) 
m=0 
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where 

mi Jq h 

Now, contributions to the contour integral in eq. (|5|) with a given power of In jhz will 
pick up the coefficients of multiple poles in the product of the expansions (0) and (§). It 
follows that going beyond leading order in (^|) introduces extra powers of M, which must 
be compensated by extra powers of 1/M in (^|), and thus extra powers of a s . Hence in 
a leading logarithm approximation we need retain only the first term in ([)]): writing the 
integrated parton distribution function Fn(/i 2 ) = f%(n 2 ) the double factorization (|5]) of 
two particle reducible contributions then becomes simply 

a N {Q 2 /ii 2 - ii 2 ) = C N {Q 2 /ii 2 - a s {n 2 ))F N {n 2 ) + 0(1/Q 2 ), (10) 

i.e. the mass factorization of the complete cross-section at high Q 2 . 

Renormalization group invariance (i.e. /i 2 independence) of the cross-section then 
implies that if 

^—Fn^ 2 ) = lN (a s (v 2 ))F N (v 2 ), (11) 
the cross-section may be written as 

a N (Q 2 /fi 2 ; fi 2 ) = Cat (1; a s (Q 2 ))T N (Q 2 , ^ 2 )F N {^ 2 ) + 0(1/Q 2 ), (12) 
where Ctv(1; a s ) = Cjy(ct s ), and the evolution factor 

rvs(Q 2 ) da 

rjv(Q 2 ,^ 2 ) = exp / ^r^lN{a), (13) 

with /3(o! s (/x 2 )) = da s ((i 2 )/d\n (i 2 . All the dependence on fx 2 may then be eliminated 
explicitly by setting F N (Q 2 ) = T N (Q 2 , fx 2 )F N (^i 2 ). 

This derivation of leading-order mass factorization may be extended to all orders in 
perturbation theory by choosing suitable mass factorization schemes (where, in particular, 
all terms in the expansion eq. (H) have the same fi 2 dependence). Eqs. (|10D- ([T3|) then still 
hold but with 

oo 

F N (Q 2 ) = f° N (Q 2 ) + C° N (a s (Q 2 ))- 1 £ C%(a s (Q 2 ))f%(Q 2 ). (14) 

m=l 

Beyond leading order this expression involves increasingly higher logarithmic moments @ 
of the unintegrated parton distribution function, so that it is in principle not possible to 
reconstruct the unintegrated parton distribution function from the integrated. 

Our main interest in this derivation is that because of the symmetry between Q 2 and 
5" 2 we may use a very similar argument to derive an energy factorization theorem valid 
at high energies, i.e. as S 2 /A 2 — > oo, rather than at high virtualities. In the high energy 
limit two particle irreducible graphs are no longer power suppressed. However they are 
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still free of logarithmic singularities [ 10 1 , at least at a fixed order in perturbation theory, 
and thus in the ultraviolet they are at most constant. Graphs which grow logarithmically 
are thus two particle reducible (in either the s or t channel), and in fact, as can be shown 
by explicit calculation, must be two gluon reducible. This means that we can put all the 
two quark reducible and two particle irreducible graphs together into a contribution to 
the cross-section which, though not perturbatively calculable, behaves as a constant in the 
high energy limit (see fig. 1). The remaining two gluon reducible contributions are then 
given by 



<J {2 \S 2 /^ 2 ^ 2 ) = — dN{^) C NM (a s (^))f NM (^), (15) 

where again the contour passes just to the right of the origin and is closed on the left, and 
the rightmost singularity is that at N = 0, the remaining singularities giving contributions 
which are power suppressed as S 2 /A 2 — > oo. 

It follows that the leading behaviour at high energies is found expanding Cnm and 
/jvm m the neighbourhood of N = 0. All the logarithmic singularities are again in the 
coefficient function: the parton distribution function is two particle irreducible and thus 
has at most a simple pole at iV = 0, so we expand it as 

oo 

f NM (n 2 )=N- 1 J2fM^ 2 )N n , (16) 



n=0 



so that 



(17) 

and in particular for n = we have simply 



/^(^) = lim/ M (y;^). (18) 

The coefficient function admits an expansion analogous to (|6|): 

oo 

C NM (a s (ii 2 )) =Nj2Cl I (a s (^ 2 ))N- n - 1 , (19) 

— oo 

so that 

oo 1 

C M (S 2 /l 2 ^ 2 ) 2? Y,n C M +1 (^ 2 ))(lnS 2 /l 2 r + 0(l 2 /S% (20) 

where C]^(a s ) begins at 0(a™) when n > 0. Because the top rung is necessarily a quark, 
which cannot generate a logarithm of S 2 , there is one more power of a s in this expansion 
than in the similar expansion (Bf). 
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Energy factorization now follows in precisely the same way as mass factorization: in 
the leading logarithm approximation we need retain only the first term in (^6|): calling this 
Fm(^ 2 ) we immediately get 

a M (S 2 /^ 2 ; a* 2 ) = C M (S 2 /» 2 ; a s ^ 2 ))F M (p 2 ) + a q M (S 2 /^ 2 ; ^ 2 ) + 0(1/ S 2 ), (21) 

where a q M contains the two gluon irreducible terms [terms (b) and (c) in fig. 1] and ap- 
proaches a constant in the high energy limit, while the first term [term(a) in fig. 1] grows 
logarithmically with S 2 . 

We may now again use renormalization group invariance of the cross-section to derive 
the longitudinal evolution equation 



» 2 -^Fm(v 2 ) = lM{a s (^))F M (fi 2 ), (22) 



where 7m (ct s (^ 2 )) is a new longitudinal anomalous dimension (as opposed to the conven- 
tional transverse anomalous dimension 7at(q: s (// 2 )) defined in (|Tl|)). The energy factoriza- 
tion (|21| ) then becomes 

a M {S 2 li?- 1?) = C M (1; a s (S 2 ))T M (S 2 , ^ 2 )F N (^ 2 ) + a q M (S 2 /^ 2 ; ^ 2 ) + 0(1/ S 2 ), (23) 

where Cm(1', a 8 ) = C^(a s ), and the evolution factor 

,(v 2 



r as(s2) da 

r M (S 2 ,ti 2 )=exp 7^7m(«). ( 24 ) 



Again these results may be extended to all orders in perturbation theory by choosing 
a suitable class of factorization schemes: the longitudinal integrated parton distribution 
function is then given by 

oo 

Fm(S 2 ) = f M (S 2 ) + CU*s(S 2 ))- 1 Y,CMMS 2 ))fUS 2 ), (25) 

n=l 

so that again the unintegrated parton distribution function cannot be reconstructed from 
it. It follows that the two integrated parton distribution functions Fn(h 2 ) and Fm(h 2 ) 
are not directly related to one another. 

Just as the usual evolution eq. ( |1~2"| ) resums logarithms of the transverse scale Q 2 , the 
evolution equation ( p3|) resums logarithms of the large longitudinal scale S 2 . Thus while 
the former evolution is asymptotically free at high virtualities, the latter is asymptotically 
free at high energies, i.e. at small x. In both cases, the resummation follows directly 
from renormalization group improvement of the dependence on the unphysical variable \j? : 
the scale of this dependence is then necessarily the large physical scale which underpins 



the appropriate factorization formula, eq. ( flOD and (|21| ) respectively. An inappropriate 
choice of renormalization scale in either case, such as for example Q 2 in a high energy 
process, results in a badly behaved perturbative approximation, signalled by a large scheme 
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dependence of the results, which is due to the fact that large logarithms are not being 
resummed. Both evolutions require nonperturbative input: in the former case a starting 
distribution Fn(Qq), and in the latter case a starting distribution Fm(Sq). However, the 
renormalization group improved factorizations ([l^) and (|23| ) imply that the final results 
are independent of the starting scales, Qq and Sq = A 2 /xq respectively. 

Before we can consider perturbative evolution in the longitudinal case we must first 
know the longitudinal anomalous dimension 7m(«s) 5 or, rewriting ( |2"2l ) as 

S 2 ^F(Q*, S 2 ) = ^1 J™ ^P(Q 2 /k 2 ; a 8 (S*))F(k 2 , S 2 ), (26) 

the splitting function P(k; a s (S 2 )) given by the inversion of the Mellin transform 

7m(« s ) = ^ r-K M P(K;a s ). (27) 
^ Jo K 

The splitting function may be calculated in the Weizsacker- Williams approximation 
by considering the emission of quasi-collinear partons in the infinite momentum frame 
(so that the transverse momenta k are much smaller than the longitudinal momenta p), 
and picking out the coefficient of whichever large logarithm we are interested in. In the 
transverse case this is related to the collinear divergence of a parton propagating in 
the t-channel with small transverse momenta k 2 , while in the longitudinal case it is the 
infrared divergence due to the emission of a soft parton of longitudinal momentum yp. 

When computing the transverse splitting function, it is sufficient to consider i-channel 
emissions, since emission in the s-channel produces no extra collinear logarithms. However, 
for the longitudinal splitting function, both s and t channel emissions can contribute: for 
purely gluonic processes the sum of the amplitudes for emission of t and s channel gluons 
(see fig. 2) is (suppressing colour factors for simplicity) 



where k = k x + ik y and the subscript ++ denotes the helicities of the emitted and final- 
state gluons: amplitudes for other helicity combinations are found by complex conjugation 
of the pertinent momenta. Evaluating the photon-parton cross section in the Weizsacker- 
Williams approximation in terms of the spin and colour averaged square of the ampli- 
tude (|28|) and the total inclusive cross section for the photon to scatter off parton C 

= ^£^£^1 (29) 



y k 2n |1 — k\ 
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where k = \k'\ 2 /\k\ 2 . The photon-parton cross section 

dcr^ thus acquires an extra energy 
logarithm due to the extra emission. The real contribution to the splitting function may 
thus be simply identified with the coefficient 2Ca/\1 — k\ of this logarithm. Virtual gluon 
emission produces a contribution to the splitting function proportional to 5(1 — k) (be- 
cause it leaves the momentum of the propagating parton unchanged), and regularises the 
singularity at k = 1, to give 



P{k) = 2C.4 



2- 1 



1 - K 



B(l - k) - P— ^- + cS(l - k) 

+ 1 — K 



(30) 



where is the step function, P indicates the principal part, + is the plus prescription 
defined in ref. ||, and c is a finite real number, which may determined to be zero by 
explicit calculation of the one loop virtual corrections, or from the knowledge of the two- 
loop transverse splitting function, as described below. 

The longitudinal anomalous dimension may then be determined from the splitting 
function (|30| ) by taking Mellin moments (|27|): 

7M (a s ) = CU— X(M) + 0(a 2 s ), (31) 

7T 



where 



\(M)= I -^(k^ + k^ -2) =2^(1)- ^(M)-^(l-M), (32) 

1 — K 



and ip(z) is the digamma function. 

Solving the evolution equation ( |2"!2| ) with the anomalous dimension ( |3"T| ) corresponds 
to iterating the s- and t-channel emissions of fig. 2, with the kinematic constraint of strong 
ordering of longitudinal momenta 1/yi << ... << l/y n « The iteration may now 
be checked by explicit computation: this is one of the advantages of using the Weizsacker- 
Williams approximation. The diagrams generated in this way are then all (cut) gluon 
ladders with rungs dresses by rainbow self-energy corrections (plus virtual corrections). 
These are the same leading logarithmic corrections as those to the unintegrated parton 
distribution function summed at fixed coupling by the BFKL equation P,|6|JT^1: x(M) 
is the Lipatov characteristic function and thus the longitudinal splitting function P(k) 
coincides with the kernel of the (forward) BFKL equation. The evolution equation (|26|), 
because it has been derived from the appropriate factorization theorem (|2~T|), performs 
instead a renormalization group resummation of these logarithms which drive the evolution 
of the integrated parton distribution function Fm(S 2 ). The perturbative approach is then 
justified self consistently because the coupling becomes asymptotically free in the high 
energy limit. 

When quarks are included, it is easy to check by a similar calculation that only the 
process of gluon emission from a quark can generate energy logarithms. Thus at the level 
of leading logarithms any initial quark must immediately turn into a gluon, in which case 
its contribution may be absorbed into the initial gluon distribution, or else be struck by 
the virtual photon, in which case it contributes to the asymptotically flat term a q M in ( |23|) . 



M 

It follows that at the leading order there is only one type of evolving parton, namely the 
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gluon. It is not difficult to see that one can choose factorization schemes in which this 
remains true to all orders. 

Finally, coefficient functions Cm(1;c*s) may be calculated at leading order by evalu- 
ating the quark box for an incoming off-shell gluon, taking Mellin transforms and setting 



N = 0. For the contributions of light quarks to F2 and Fl this gives [11,1 



2 « s 2(2 + 3M-3M 2 )r 3 (l-M)r 3 (l + M) 

M{ ' s) 2tt R M(3-2M) r(2 + 2M)r(2-2M) + { sJ ' 

n • „ ^ - ^£r 4(1 -M) r 3 (i-M)r 3 (i + M) , (33) 

M[ ' sJ " 2tt fl 3-2M T(2 + 2M)r(2-2M) + 
coefficient functions for heavy quark production may be deduced from the results in 



ref. HTg,^7|. Whereas at leading order the M-dependence of any particular coefficient 
function can be absorbed in a redefinition of the starting gluon distribution Gm(S 2 ), and 
has thus no physical meaning, ratios of coefficient functions are physically meaningful. 

The complementary role of the longitudinal and transverse renormalization group 
equations fl22|) and (jll|) may be further elucidated noting that they are connected by a 
duality relation. If for simplicity the coupling is held fixed and the equations are truncated 
at leading order, they take the form 

' ;*m(S 2 ) = a sX (M)a M (S 2 ), t^^<?n{Q 2 ) = a g( p(N)a N (Q 2 ), (34) 



where the leading order pure gluon transverse anomalous dimension 7^(ct s ) = a s (j)(N), 
with a s = — a s , while the corresponding longitudinal anomalous dimension is given by 
(|31~D. Taking a second Mellin transform of both equations, we then find that 

Na NM = o, s x{M)a NM , Ma NM = a s <p(N)a NM , (35) 



respectively. So both equations can be viewed as defining a trajectory in the N-M plane, 
each of which approximates at leading order in a s the "true" trajectory (as given by a 
hypothetical all-order computation), the former when N is small, the latter when M is 
small. Both equations will only be valid together if both iV and M are small: since both 
X and (f) have a simple pole with unit residue at the origin, both equations then reduce 
to NM = a s , which corresponds (at fixed coupling) to the double logarithmic or double 
scaling approximation |T8| , |i~9|| . 

This however also implies that we may turn the longitudinal renormalization group 
equation into a transverse renormalization group equation, and vice versa, by defining 
leading singularity anomalous dimensions implicitly as 

1 = fx(7^(f)) ) l = f^(7f S (f)), (36) 
in terms of which the two evolution equations (|35|) may be written as 

Ma NM = lT S (if)°NM, No NM = ji S (jj)°NM, (37) 
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respectively. Thus the leading singularities in N can also be summed up explicitly 
using a transverse renormalization group equation, provided the anomalous dimension 
7jv(os) = 1t S (jt^ * s use d 7 which which adds up the leading singularities in N in the 
transverse anomalous dimension to all orders |Il, pO|,21[] ; conversely the M singularities 



may be summed up using a longitudinal renormalization group equation. 

This duality can thus be used to provide powerful cross-checks on the transverse 
and longitudinal anomalous dimensions, since knowledge of the leading order transverse 
anomalous dimension determines the leading singularities of the longitudinal anomalous 
dimension to all orders, and conversely. For example the result c = in (|30|) is equivalent to 
the vanishing of the leading singularity at NLO (i.e. the 0{a 2 s /N 2 ) term in the transverse 
anomalous dimension ||12| ). It is now also clear that the failure of attempts |2!| to extend 
the reach of the transverse evolution equations to high energies (JTj) by simply adding the 



leading singularities 7^ s (ff ) to a fixed order truncation of 77v(a s ) is due to the fact that 
while the mass singularities are being resummed (because the coupling runs with Q 2 ), the 
energy singularities are merely summed f21 ]. This explains the large scheme dependence 



in such calculations |23|1 , and ultimately their failure to agree with experiment |7|,|24|]. 

We can now discuss some of the consequences of our results for the behaviour of 
cross-sections at asymptotically high energies. The structure function at leading order is 
determined by solving the evolution equation (^3j) with the anomalous dimension (^T[), and 
combining with the appropriate coefficient function (|33|) : 

pe+ioo ,rp-\M /In <? 2 / A 2 \ -y 2 Y(M) 

F2(x ' Q2) = 2^i . dM iw ^ e2) ^ (1;as( ^ 2) Khrl^) g °m+ f 2(q% 

° (38) 

where (e 2 ) is the average of the light quark charges, 7 = ^ACa/Po, and /3q is the leading 
coefficient of the /3-function. The a priori unknown inputs at the starting scale Sq are the 
initial integrated gluon distribution G°(Q 2 ) (the only evolving longitudinal parton distri- 
bution) and the unevolving quark piece F^Q 2 ). To determine the asymptotic behaviour 
at high energy we may take a dipole expression Q 2 /(Q 2 + A 2 ,) for G°(Q 2 ), which has the 
correct behaviour at both large and small Q 2 . Its Mellin transform is then 

G O f .JW£!) M (39) 
sin7rM VAq/ 

At high energies (and thus large rf) the integral (|38J) will be dominated by the contri- 
bution from the saddle point at M s , where 

t + 7 V(M s )7? = 0, (40) 

where t = lnQ 2 /Ao (note that the nonperturbative input sets the scale for Q 2 ) while 

In / fV^ 1 

7] = In lnS y A 2 = l n hf^- Since x'(M) increases monotonically as M increases from zero 
to one, the precise position of the saddle point will depend on the ratio t/rj. 
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Consider first the high energy limit rj — > oo, t/rj — > (when t is fixed this is the 
so-called Regge limit - see fig. 3). Then expanding x'(M) about M = i, we have M s = 
\ — t/28(37 2 ?7 + • • •) whence the asymptotic behaviour 

F 2 (x,Q 2 ) ~ Q 7Vof^J 1,-5 flu-) + (41) 

valid in the region rj ^> |t|/28(37 2 . Subleading corrections to this behaviour are down by 
powers of \t \/rj, while corrections due to higher orders in perturbation theory are suppressed 
by powers of In 1/x. 

It follows that in the high energy limit F 2 grows logarithmically at a universal rate, 
quite independently of the starting conditions, since eventually the first term in (|4l| ) must 
always overcome the second. Furthermore, in this limit the Q 2 dependence is also universal. 
The universality is a direct consequence of the fact that their is only one evolving parton at 
high energies. This parton can be viewed as playing the role of the pomeron: it generates 
a universal rise in the high energy limit. However the rise is only logarithmic, rather than 
powerlike, which is compatible with the logarithmic nature of unitarity bounds [§. 

Several asymptotic predictions can be derived from this universality, and its applica- 
bility to other structure functions. For example using (|3^) we see that asymptotically F 2 
and Fl grow in the same way, and in fact 

F L (x,Q 2 ) _ 2_ 
*™F 2 (x,Q 2 ) " 11- {42) 

More generally we can derive sum rules 

I 00 dQ2 (0^~ M F (r O 2 ) - 2 + 3M(l-M) r dQ 2 2 M 2 
J -W iQ) F2(X ' Q) " 2M(1-M) Jo ~^ {Q) Fl ^ Q ^ < 43 > 

valid for < M < 1. Similar sum rules may be derived for other inclusive structure 
functions. Corrections are now down by powers of a s (S 2 ), rather than powers of rj. 

Let us now turn to the opposite limit, i.e. that in which as r\ — > oo, \t\/rj — > oo. In 
the high virtuality limit t/rj — > +oo the position of the saddle point shifts to just above 

the origin: M s = ^ — 2^3^ + . . ., where we write a = \/\t\r) and p = \f\t\/r), (similarly 

to [jT^|, but with the roles of - and Q 2 interchanged). We then find that asymptotically 

F 2 (x,Q 2 ) ~ Af+a-^pe^-^+F^Q 2 ), (44) 
1 

valid whenever p 3> (2^3)37. This is even softer than the universal high energy behaviour 

(PI). Similarly in the low virtuality limit t/rj — > —00, i.e. the photoproduction limit 

4 

(see fig. 3), the saddle point shifts to just below unity: M s = 1 — ^ + 2(3^ + . . ., and 
asymptotically 

F 2 (x 1 Q 2 ) ~ M-(Q 2 /K 2 )a- l 2 P 2 e 2 ^-^ + F|(Q 2 ). (45) 

<3 2 ^o 
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The extra factor of Q 2 now ensures the finiteness of the perturbative contribution to the 
photo-production cross-section, while the growth with energy is softer than any power of 
lnS' 2 , and thus inside the Froissart bound 0, as it should be. For Fj_, an extra power of 
Q 2 is generated by the coefficient function provided that the quark mass is kept finite. 



Both the results (|4J) and fl45| ) must however be treated with caution, since whenever \t\ 
becomes large we start to leave the the region of validity of the longitudinal renormalization 
group equation: logarithms of Q 2 become as important as logarithms of S 2 . This can 
be seen in the perturbative expansion of the anomalous dimension 7m(o s ): as M — > 
higher order terms behave as (a s /M) m , generating large ultraviolet logarithms, while as 
M — > 1 they behave as (a s /(l — M)) m , generating large infrared logarithms, and in both 
cases leading to the breakdown of the perturbation series. This is precisely analogous to 
the breakdown of the transverse renormalization group equation as x — > and x — ► 1 
respectively. 

The contributions from the infrared regions (i.e. the Sudakov region and the pho- 
toproduction region) are in both cases suppressed by the fact that they are kinematic 
boundaries, and thus the parton distribution functions are small. The ultraviolet regions 
instead are complementary, since the behaviour (]4^) matches smoothly into the double 
scaling behaviour [1^,|19| . In particular the asymptotic behaviours (f|l]) and (|4[) provide a 
soft boundary condition to the transverse renormalization group equation at small x (i.e. 
rising more slowly than any power of 1/x), so as Q 2 increases the double scaling rise even- 
tually dominates. This explains why double scaling is so successful phenomenologically 

In conclusion, we have discussed a new factorization theorem, valid at asymptotically 
high energies but moderate values of Q 2 , which correctly resums all high energy logarithms. 
We found that there is then a consistent partonic picture of high energy scattering in which 
the partons are asymptotically free. There is only one type of evolving parton, which 
may be thought of as a gluon. As it evolves it generates a universal growth of structure 
functions at small x reminiscent of that traditionally described by the pomeron, but much 
softer (logarithmic), and thus consistent with known unitarity limits. This might have 
far reaching implications not only for our understanding of other inelastic lepton-hadron 
processes (such as diffractive and rapidity gap events), but also for photon- photon and 
purely hadronic processes. 
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Figure 1: Contributions to the forward amplitude: a) two gluon reducible, b) two gluon 
irreducible but two quark reducible c) two particle irreducible. At high Q 2 a) and b) are 
leading twist and evolve, while c) is higher twist (power suppressed). At high S 2 none 
of them are power suppressed, but only a) evolves, while b) and c) are asymptotically 
constant. 
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Figure 2: Diagrams for (a) t-channel and (b) s-channel gluon emission contributing to 
the splitting function. For each line the parton's energy (in units of p) and transverse 
momentum are indicated. 
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Figure 3: Regions and limits in the (Q 2 , ^) plane. The shaded area denotes the region 
of nonperturbative dynamics where there is no large scale. 
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